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We study an ultra-cold gas of A'^ bosons trapped in a one dimensional Af-site optical lattice 
perturbed by a spatially dependent potential g ■ x-' , where the unknown coupling strength g is to 



be estimated. We find that the measurement uncertainty is bounded by Ag oc . For a 

typical case of a linear potential, the sensitivity improves as M^^, which is a result of multiple 
interferences between the sites - an advantage of multi-path interferometers over the two-mode 
setups. Next, we calculate the estimation sensitivity for a specific measurement where, after the 
action of the potential, the particles are released from the lattice and form an interference pattern. 
If the parameter is estimated by a least-square fit of the average density to the interference pattern, 
the sensitivity still scales like for linear potentials and can be further improved by preparing 

a properly correlated initial state in the lattice. 
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I. INTRODUCTION 

As scientists gain insight into the "nano" scale in both 
the space and time domain, it becomes increasingly rele- 
vant to develop ultra-precise measurement devices. Very 
prominent among those are interferometers, where the 
parameter to be measured is mapped into a phase differ- 
ence 9 (or optical path length difference) between two (or 
more) arms of the device. By recombining the outputs of 
each arm, an interference figure appears, which contains 
information about the parameter to be estimated. 

The key quantity which describes the performance of 
an interferometer is its precision A9, often referred to 
as the sensitivity. In a typical interferometer, where the 
signal which is used to determine 6 comes from a collec- 
tion of N independent particles, the sensitivity scales at 
the shot-noise level, i.e. A6 cx -i — i=, where m denotes 
the number of independent repetitions of the measure- 
ment. This expression tells that the sensitivity can be 
improved by either increasing the number of trials m or 
the number of probe particles N. A third way of de- 
creasing Ad is by replacing the uncorrelated ensemble of 
atoms with a usefully entangled state. Non-classical cor- 
relations can be employed to beat the Shot-Noise Limit 
(SNL) and ultimately reach the Heisenberg Level (HL), 
where A9 cx — ^ 4? , which for large N can be a significant 
improvement over the SNL 

A paradigmatic example of an interferometer which 
benefits from the use of entangled states is a Mach- 
Zehnder Interferometer (MZI), where the signal coming 
from two arms at the input is first recombined/split using 
a beam- splitter, then the phase difference is imprinted 
between the two arms, and finally another beam-splitter 
maps the phase difference into the intensities in the out- 
put ports. When the MZI is fed with an entangled spin- 
squeezed state and the phase 6 is estimated from the 



measurement of the population imbalance between these 
ports, the Sub Shot- Noise (SSN) sensitivity AO < 

can be acheived. 

While originally implemented with photons, in recent 
years interferometers fed by atoms have attracted a lot 
of interest and will be the object of the present work. 
These interferometers proved promising for the precise 
determination of electromagnetic [l-Q and gravitational 
interactions. Moreover, their performance can be 
improved with the recentl y c reated non-classical states 
of of matter (9l-[T6|. In [Tsi - llH a spin-squeezed state was 
generated by reducing fluctuations between the sites of 
a double-well potential, while in [12i, .16 18] internal de- 
grees of freedom were employed. In the experiments men- 
tioned above, a two-path (or two mode) usefully entan- 
gled quantum state was prepared. In this work, we con- 
sider instead a multi-path atom interferometer, involving 
more than two internal/external degrees of freedom. 

The idea that the growing number of modes could im- 
prove the performance of an interferometer has been first 
discussed in the context of optical devices [l^, [l^] . They 
involve the multi-mode version of a beam-splitters, and 
allow in general to measure a vector phase [2l|. In the 
particular case where the phase difference is the same be- 
tween each pair of "neighboring" paths, different meth- 
ods have been used to derive bounds for the phase sensi- 
tivity. In (20| . for an input coherent state, a momentum 
formula approximation allowed to predict a sensitivity 
AO cx M~^, where M is the number of ports. This is a 
result of a multiple interference between different paths, 
and does not violate the HL "25] . The gain is also robust 
against losses (20j , differently from the improvement due 
to entanglement present in the input state. The scaling of 
the sensitivity with the number of modes and particles 
has also been heuristically derived for the specific case 
where coherent or number states are used in input [23| . 
Recently, the sensitivity for rotation measurements using 
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a ring interferometer with ultracold atoms inside a lattice 
has been studied in [24] . Numerical determination of the 
ground state allowed to calculate the Quantum Fisher 
Information (QFI) for specific states, which, as we will 
discuss in the following, provides the ultimate sensitiv- 
ity bound, independently of the measurement and esti- 
mation strategy employed. On the experimental side, a 
multimode atom interferometer was implemented in [25| , 
where a single Bose-Einstein Condensate was coherently 
split among five Zeeman sublevels using a radio frequency 
pulse. The relative phase between the atoms occupy- 
ing different substates was imprinted with an external 
magnetic field. Finally, the cloud was recombined with 
another pulse and an increased interferometric signal - 
with respect to the two-mode case - was observed. A 
multi-mode nonclassical state of ultracold bosons inside 
an optical lattice has also recently been created [26j 

In this manuscript, we consider a multi-path interfer- 
ometer realized with a dilute gas of ultracold bosons 
trapped inside a one-dimensional optical lattice [27l - l29| . 
Each of the M lattice sites plays the role of a possible 
interferometric path, whose optical length can be modi- 
fied by adding a spatially dependent force which would 
be then the object of our precision measurement. More 
precisely, for a potential V{x) = gx^ , we derive the sensi- 
tivity for the estimation of the parameter g. First, using 
the Quantum Fisher Information (QFI), we derive the 
best possible precision to be Ag oc ^/^^jy^^.p ,!-) j indepen- 
dently of the quantum state, output measurement and 
estimation strategy. Then, choosing a specific class of 
states which are symmetric with respect to exchange of 
lattice sites, we obtain the sensitivity Ag oc , — 

for M ^ 1. Here, f{N) is a two-mode analog of the 
enanglement measure, which can vary from TV for the 
superfluid state of the lattice up to N"^ for a strongly en- 
tangled state. Note that when the potential V{x) acts on 
the system, one can naturally define a phase — gx^t/h, 
which is related to the coupling constant, the duration 
of the interaction t and the lattice spacing xg. In spirit 
of the two-mode interferometry, one can thus speak of 
the phase estimation, rather then determination of the 
coupling constant g. Throughout this work, we use both 
Ag and A0, and keep in mind they are simply related by 
Ag = (n/te^)A6l. 

We finally consider a specific measurement and esti- 
mation strategy which could be applied in most of the 
current ultracold atom experimental setups: after un- 
dergoing the action of the force, the atoms are released 
from the lattice, form an interference pattern, and the 
parameter is extracted by a least-squares fit of the aver- 
age density. Using rigorous results from the estimation 
theory |30i] , we derive the sensitivity for this interferomet- 
ric protocol, and show that it scales with M in a same 
manner as in the case of the optimal measurement. The 
use of the spatial interference pattern, obtained through 
the expansion of the cloud, substitutes the beam splitter. 
The in-situ implementation of the latter in an optical lat- 



tice, in analogy to the double- well implementation, seems 
a much more demanding task. 

The manuscript is organized as follows. In Section 
im we formulate the framework for analyzing the perfor- 
mance of the multi-mode interferometers. In Section [III[ 
we employ the notion of the QFI to provide some ulti- 
mate bounds for the multi-mode interferometry. First, in 
Section UlIAl we identify the states, which allow to reach 
the best possible scaling of /S.9 with the number of sites 
M and linear potentials V{x) — gx. Then in Section 
IIIIBl using the QFI, we show how Ad scales with M for 
states which are symmetric with respect to interchange of 
any two sites. In Section UlI CI we generalize these results 
for non-linear potentials V{x) = gx^ . Finally, in Section 
IIVI we show how the scaling with the number of sites 
can be achieved for a particular detection scheme and 
identify the states usefully entangled for the multi-mode 
interferometry. Some details of analytical calculations 
are presented in the Appendices. 

II. FORMULATION OF THE MODEL 

We begin the analysis of multi-mode interferometry 
by introducing the quantum state Itp) oi N bosons dis- 
tributed among M modes. The most general expression 
for \tP) is 

N N-711 Af-ni-...-riM -2 

1'^) = E E ■•• E (la) 

m— 712—0 UM-l—O 

|ni,n2, . . .,nM-i,N- ni - . . . - um-i), (lb) 

where the sum runs over occupations of M — 1 wells. 
Since the total number of atoms is set to be A^, the pop- 
ulation of the M-th mode is determined by other M — 1 
occupancies. The coefficient Cm. ..um-i is the probabil- 
ity amplitude for having ni . . . hm-i atoms in the M — I 
modes and A^ — (ni -f . . . um-i) in the last one. 

For future purposes we introduce a shortened nota- 
tion, where the set of M — I indeces (ni . . . riM-i) is 
represented by a vector n, and the state ([T]) is written as 

\,p)=J2Cn\n). (2) 

n 

The complete description of quantum properties of the 
system includes the M-mode field operator which reads 

M 

•^{x)=J2M^)ak. (3) 

k=l 

Here, Ofc is a bosonic operator which annihilates an atom 
occupying the mode-function ^pk{x) localized in the fc-th 
well. 

Let us for now restrict ourselves to the case where 
a constant force, corresponding to the linear potential 
V{x) = gx, is added to the lattice potential (the more 
general case will be consider in Section IIII Cp . Here g is 
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the couphng constant - the parameter which we want to 
estimate. 

The Hamiltonian corresponding to the added poten- 
tial, in the second quantization, reads 



III. SENSITIVITY FOR THE OPTIMAL 
MEASUREMENT 

A. Ultimate scaling 



M 



H ^ ey^ k ■ fik, 



(4) 



fc=i 



where fik — a\ak is the on-site atom number operator 
and e = g ■ xq, where xq is the lattice spacing. 

The external force interacts with the trapped gas for a 
time t, which transforms the field operator ([3]) as follows 



^(a;|6l) 



iOh 



-ike 



(5) 



fc=i 



where 9 — et/h and h is the generator of the phase-shift 
transformation 



M 

h = y^ k ■ fik- 



(6) 



fc=i 



Equation ([S|) shows that the perturbation imprints a 
phase between two neighboring sites, while the phase 
between the two most distant wells is {M — 1)9. 

The interferometric sequence consists thus of a simple 
phase shift acting on the input state ([2]) and could be 
implemented as follows. First, is the preparation of the 
input state ([2]) in an unperturbed optical lattice (g — 0). 
Then the perturbing potential is suddenly turned on on a 
time scale shorter than the lattice hopping time, but slow 
enough to avoid creating too many excitations. In this 
way, the state does not follow the change adiabatically, 
and the evolution operator is exp{—ih9). The goal of this 
work is to derive the sensitivity of the estimation of or, 
equivalently, of the coupling constant Ag = {h/txo)A9. 

Given the interferometric operation defined by 
Eqs. (O,®, the sensitivity A9 still depends on: a) the 
choice of the input state lip), b) the choice of the measure- 
ment to be performed on the final state exp{—ih9)\'ip), 
and c) the choice of the estimator, that is, how to in- 
fer the value of the parameter given the measurement 
outputs. In Section [ml we will derive the optimal sensi- 
tivity, that is, assuming that both the measurement and 
the estimator are the best achievable. This can be done 
using the notion of the QFI. The latter allows to calcu- 
late the optimal sensitivity for a given input state. In 
Sections IIII Al IIII C[ by further optimizing the QFI over 
all possible input states, we derive the ultimate scaling 
of the sensitivity. In Section IIII Bl instead, we restrict to 
a particular family of input states which are symmetric 
with respect to the exchange of two lattice sites. Finally, 
in Section IIVI we derive the sensitivity for a particular 
measurement and estimator, using the familty of input 
states used in Section IlII Bl 



As anticipated, using the notion of the QFI, denoted 
here as Fq, we can calculate the sensitivity for the best 
possible measurement and estimator. Once we know Fq, 
we can indeed employ the Cramer-Rao Lower Bound 
(CRLB)[3l|: 



A'^9 > 



1 1 

m Fn 



(7) 



where m is the number of independent repetitions of 
the measurement used to estimate 9. The right side of 
Eq. ([7]) is the sensitivity optimized over all possible mea- 
surements and estimators. 

The QFI yet depends on the interferometric operation 
and the input state. For pure states, the value of Fq 
is given by four times the variance of the generator of 
the interferometric operation, in our case the operator 
h generating the phase shift, calculated with the input 
state of the system and reads 



= 4A'^h. 



(8) 



The QFI can be maximized with respect to IV^) and the 
upper bound for Fq is provided by a simple relation 



A^h < i(A„ 



X • ^2 



(9) 



where Amin/max are the smallest /largest eigenvalues of 
the operator h. The two eigenvectors of h with extreme 
eigenvalues are 

I V'min) = |iV, 0, . . . , 0) ^ h\^mXu) = V-min) , (10a) 
IV-max) = |0, 0, . . . , iV) ^ /l|7/;,„ax> = MA^I^-max)- (10b) 

By combining Equations ([9]) and (ITOl) we get Fq < 
N^{M — 1)2. The inequality is saturated for the 
NOON state ^ (IV'min) + IV'max)) for which, according 
to Eq. d?]), the sensitivity is bounded by 

1 1 1 



A^9 > 



mN^ {M - 1)2 



(11) 



This expression is the first important result: when the 
perturbation is a linear potential (constant force), the 
best sensitivity reachable scales like - a kind of 

"Heisenberg scaling" with the number of wells. In section 
IIII C[ we will derive the equivalent bound for nonlinear 
potentials. 

As anticipated in the introduction, this scaling is a 
purely geometrical - and thus not quantum - effect. 
While the Heisenberg scaling with the number of atoms 
N originates from a macroscopic quantum superposition 
present in the NOON state, the scaling with M is a re- 
sult of an accumulated phase shift between two extreme 
sites, as discussed below Eq. ©. 
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B. Scaling for symmetric states 

Our next step is to remove the optimization of the QFI 
over ah possible input states ji/'), and concentrate on a 
particular family of states which are more relevant for 
current experimental setups with ultracold atoms. We 
assume thus that N atoms are distributed among M sites 
symmetrically, i.e. the coefficient Cn is invariant upon 
exchange of any two indices. A prominent member of 
such family of states is the superfluid state [131, where 
each atom is in a coherent superposition, i.e. 



|0). 



(12) 



This state can be written in the form ^ with coefficients 



VM^ ^ni\...nM-i\{N -m 



■ ~ "-M-l)! 

(13) 

which clearly posses the aforementioned symmetry. An- 
other example is a symmetric Mott insulator state [27| . 
where each site is occupied by N /M atoms. Below we 
calculate the sensitivity using a generic symmetric 
state without taking any specific value of C^'s. 

To this end, we calculate the QFI using Eq. ([S]) and the 
phase-transformation generator h for a linear perturbing 
potential from Eq. The average is equal to 



M 



k=l 



1)N, 



(14) 



where in the last step, due to the imposed symmetry, we 
used (fifc) = jj. Calculation of the average of involves 
a few intermediate steps, which are presented in detail 
in Appendix [Xj The final result is that, for symmetric 
states and a linear perturbing potential, the QFI reads 



(15) 



where A^n = (fi^) — (n)^ is the on-site variance of the 
atom number. The value of A^h can be easily calcu- 
lated for some particular states. For instance, the un- 
correlated super-fluid state ([T2|) gives A^n — ^^N, 
and as a consequence the QFI scales at the SNL, i.e. 
Fq = i(M^ — 1) • A^. A strongly entangled symmetric 
NOON state 



IV': 



NOON/ 



|iV,0,...,0)-f |0,7V,...,0)-f ... + |0, 



,N) 



M 



(16) 

gives A^n = ^^jj^-N^, which gives the Heisenberg scaling 
of the QFI, Fq = i(Af2 - i) . TV^. 

These two examples suggest that for all symmetric 
states the on-site variance has a particularly simple form 
A^n = ^^/(iV), where f{N) does not depend on M. 
Indeed, the above form of A^n can be explained as fol- 
lows. When the "volume" of the system ~ in this case 



the number of sites M - tends to inflnity, the atom num- 
ber fluctuations in one site tend to zero. On the other 
hand, when the whole system consists of only one site, 
the fluctuations are zero, since the number of atoms in 
the system is flxed. The only possible dependence of the 
variance on M, which behaves correctly in both these 
regimes and also accounts for the scaling of the average 
value, (n)2 = is thus A^n = I^f[N). Putting 
this into Eq. ^ gives Fq = - 1)/(A^) and the 

sensitivity 



A^ 



> 



1 1 



mf{N) AP -1' 



(17) 



The function f{N) can be interpreted as a measure of 
the entanglement present in the system, related to the 
increased on-site atom number fluctuations. In the two- 
well case M = 2, such usefully entangled states which 
give A^n> N are called "phase-squeezed" [sl], because 
the growth of the number fluctuations is accompanied by 
a decreased spread of the conjugate variable, which is the 
relative phase between the wells. 

Most important, the scaling with M of the sensitivity 
for symmetric states (jl7|) is only slightly worse than the 
ultimate scaling (fTTj) obtained for the best possible state, 
and becomes the same for large M . 



C. Scaling for non-linear potentials 



In this part, we generalize the results of section UlI Al 
for non-linear potentials V{x) = gx\ where j € M. We 
assume that the first lattice site is placed at x — and 
define the interferometric phase as 9 = igx^/h. Accord- 
ingly, the phase imprinted on the k-th site is kW. From 
this, we can easily construct the generator as in Eq. ([6]) 
which is equal to 



rik- 



fc=i 



The QFI can be again bounded by Eq. ([9]), giving 



A^e > 



m (Ai, 



^min J 



(18) 



(19) 



The extreme eigenvalues for j > are equal to Amin = N'^ 
for |-0min) and A^ax = A-P for IV'max) while for j < 
the minimal and maximal eigen-values are exchanged. 
The states IV'min) and ji/'max) are still the ones defined in 
Eq. PH)) . The best possible sensitivity is thus 



A^0 > 



1 



1 



1 



m {Mi - 1)2 N^' 



(20) 



The behavior with M of the r.h.s. of the above equa- 
tion for various j is show in Fig.[TJ When j > 0, A^O con- 
tinuously improves with increasing M , and the larger j is, 
the quicker is the improvement. This would be the case - 
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m-N-A-'e 

4H r 




FIG. 1. Ultimate scaling of the sensitivity for various 
non-linear potentials given by Eq. (|20p . Red squares, blue 
circles, black diamonds and green triangles correspond to 
j = —2,— 1,1,2, respectively. The black dashed line denotes 
the sensitivity of a standard = 1, M = 2 two-mode interfer- 
ometer. 



for instance - for the quadratic potential, i.e. j = 2. On 
the other hand, when the potential decays with distance 
j < 0, as is the case of the Casimir-Polder-type interac- 
tions which can exhibit the j = — 3 or j = —4 scaling 
[33I Is^ l , there is no substantial gain from the multi-mode 
structure of the interferometer. This happens because 
already the phase imprinted on the second well, which is 
equal to 2~^^^6, is negligible when compared to 6. In this 
case A^6 > and by increasing M the sensitivity 

saturates to the value corresponding to a linear potential 
and two wells, since effectively only the first well experi- 
ences the action of the perturbing potential. 

Apart from the ultimate scaling, it is interesting to de- 
rive an expression for the sensitivity for symmetric states, 
as was done in Section IIII Bl for linear potentials. The 
derivation does not change substantially with respect to 
the line of reasoning presented in Appendix [Aj with the 
only difference that the summations over the sites in- 
volve instead of k. The resulting sensitivity bound for 
symmetric states reads 



(r) 

where Hn 



1 1 



M2 



fiN) 



(21) 



^IJ^j^ k^'' is known as the Harmonic num- 
ber. For j > and a large number of lattice sites, the 
sensitivity bound can be approximated by 



> 



1 1 



1 (l-fj)2(i + 2j) 



m f{N) 



(22) 



which shows that the improvement on A9 is pre- 

served for symmetric states. 



We underline that - via Eq. ([20)) - the multi-mode 
interferometer could serve as a sensitive tool to determine 
the spatial dependence of the perturbing potential, which 
is a relevant task for instance in the determination of the 
Casimir-Polder force or deviations from the Newton law 
of gravity at the micrometer scale. 



IV. SENSITIVITY FOR THE FIT TO THE 
INTERFERENCE PATTERN 

The QFI was used above to set the lower limit for 
the sensitivity by optimizing of A^9 with respect to the 
output measurement and estimator. Below we demon- 
strate that the scaling of the sensitivity with M, which 
is present in Eq. (jl7p . is preserved even for a particularly 
simple phase estimation protocol, although the sensitiv- 
ity does not saturate the QFI. Also the scaling with N is 
preserved for the superfluid and phase-squeezed states. 

We assume that after the system acquires the phase 
as in Eq. ([5]), the optical lattice is switched off and the 
wave-packets ipkix) - which were initially localized at 
each site - expand and in the far field form an interference 
pattern. Positions of atoms are recorded and the phase is 
estimated by fitting the density p(a;| 6*) = {i'Hx\9)i'{x\0)) 
to the acquired data. In our past work [30i] we have 
shown that when this estimation protocol is employed, 
the sensitivity is given by 



^'9 = 



1 fi + C 
m ' 



(23) 



where Fi is the one-body Fisher information which reads 



^^1 



dx 



1 f dp{x\9) 
p{x\9) V d9 



(24) 



The coefficient C depends both on the one- and two-body 
correlations, 

fdx fdyG^^\x,y\9)dg\og[p{x\9)] de log [p{y\9)] , 



(25) 

where G^^\x,y\9) = {x\9)i''< {y\9)^ {y\9)i' {x\9)) . 
Our goal is to analyze how the sensitivity (1^ depends 
on both the number of sites M and the non-classical cor- 
relations present in the state 

In section IIII Bl we have seen that the usefully entan- 
gled states for the present interferometer have enhanced 
number fluctuations. In particular, the NOON state pro- 
vides the best scaling with the number of particles N. 
This statement, however, assumes that we are always 
able to perform the optimal measurement, which in gen- 
eral depends itself on the input state. In this section in- 
stead we fix the measurement (as described above), and 
study the behavior of the sensitivity for different states 
characterized by enhanced number fluctuations (we will 
see that such states are still the useful ones also for the 
particular measurement considered) . For two lattice sites 
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(M = 2), one can easily generate a family of such use- 
fully entangled states by numerically finding the ground 
state of the two- well Bose-Hubbard Hamiltonian 



BH 



(26) 



with U < 0. Here, Ej is the Josephson energy, related 
to the frequency of the inter-well tunneling, and U is the 
amplitude of the on-site two-body interactions. This pro- 
cedure for generating the input states is also relevant for 
the setup we consider, since it will correspond to adia- 
batical preparation in the lattice. To perform an analo- 
gous calculation for higher M is very demanding, since 
the dimensionality of the Hilbert space grows rapidly, as 
indicated by multiple sums present in Eq. ([1]). This dif- 
ficulty can be overcome by replacing the true state ([T]) 
with 



M oo 
k—l rtk—O 



(27) 



This is the Gutzwiller app roximation of the quantum 
state of the bosons (see [211 ^^.d references therein). The 
number of atoms is not fixed anymore, and thus each site 
is independent from the others. Although such state does 
not recover correlations among the wells, it can be used to 
mimic non-classical correlations by increasing the on-site 
atom number variance A^n. In analogy to the fixed-A^ 
case ([2]) we assume that the average number of atoms per 
site is (n/c) = so that the average number of atoms in 
the whole system is N. 

We express the field operator 'i>{x\9) in the far field. 



M 
k=l 



(28) 



We assumed that initially all the wave-packets V'fc (x) had 
identical shape but were localized around separated sites 
of the optical lattice. The function ip is the Fourier 
transform of ipk (x) , which is common for each well, and 

(fi is the atomic mass and t is the expansion 



time). The phase factor is defined as 



(29) 



where xq is the lattice spacing first introduced below 
Eq. Q. 

The density p{x\d) obtained using Equations (l?fl) and 
is 



The density is put into Eq. (l24l) . giving 



Fi=N dx 



' 1 [ df{x\9) 
fix\0) y 89 



(32) 



When the density of the interference pattern consists of 
many fringes, the envelope \ip\'^ varies much slower in x 
then the function f{x\9). An integral of a slowly varying 
envelope with a quickly oscillating periodic function can 
be split into a product of an integral of the envelope alone 
times the average value of the quickly changing f{x\9). 
Since the envelope is normalized, we immediately get 



2-K /(0) \ d<\> ) 



(33) 



where was defined in Eq. (|29|) . This integral depends 
on the quantum state through the average value (a), as 
can be seen from Eq. pOI) . Apart from the analytical 
case Af = 2, its value has to be calculated numerically. 

Now we have to calculate the coefficient C, which in 
turn requires to compute the second-order correlation 
function. After a lengthy calculation we obtain (see Ap- 
pendix |B] for details) 



C 



~2\2 



-2^(a)^-2(a^)(a)2 



h (34a) 



(34b) 



where the two integrals are given by 

h = g f ^ a, log [/(0)] de [fm (35a) 

(35b) 



cos[A'f ((/) -I- (j)')] 



1 



cosh 



and 



1-1 



^ a,log[/(0)]ae[/(0')] (36a) 



sin [f f 


sin [f 




sin [ 


f(0 + '^')] 


sin 


2 


sin 


0' 
2 


sin 


"0 + 0'' 

2 





(36b) 



As a first test of Equations (|33p and ([M]) . in the next 
section we show that for M = 2 we recover a known 
analytical expression for the sensitivity obtained using a 
two- mode state with fixed N ISOll. 



p{x\9) = {^A^\x\9)i>{x\9m = N\ij (^) p • f{x\e), 

(30) 

where jV^p is the envelope of the interference pattern, 
while the fringes are contained in 



f{x\e) = i + 



M 



(31) 



A. Fit sensitivity with two lattice sites 

To make a direct comparison, we start by recalling the 
main result of our previous work (3Q| . where we focused 
on phase estimation with a fixed- A'^ two-mode state \4>) = 
J2n=o Cn\n,N — n). In the resulting A'' qbits space, every 
possible interferometric operation can be mapped into a 
rotation, with the angular momentum operators Jx = 
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i(at6 + ab'f), Jy = Jj(at6 - afot) and = i(aU - 
We have shown that when the phase is deduced from the 
fit of the one-body density to the interference pattern 
formed by two expanding wave-packets, the sensitivity 
P3p can be expressed in terms of two quantities. One is 
the spin-squeezing parameter [s^ , which is a measure of 
entanglement related to the phase-squeezing of and 



reads ^0 = (j^y 
interference fringes, v 



A'9 = 



The other is the visibility of the 



jf{Jx)- The sensitivity 



mN 



(37) 



on-site atom number fluctuations. In both cases, the sen- 
sitivity has an optimal point, where the entanglement is 
balanced by the fringe visibility. Clearly, the Gutzwiller 
approximation gets worse when the correlations between 
different wells increase (that is, going further to the left 
in Fig. (2), but for the usefully states we consider and for 
this particular choice of the measurement, it is exact for 
all practical purposes. 

Having proven its usefulness, in the following section 
we employ the Gutzwiller ansatz (P7)) plus the Gaussian 
approximation on the coefficient of the product state (j39p 
in order to analyze how the sensitivity ((23)) behaves for 
higher M, where results from a fixed- model are un- 
available. 



results from the interplay between the growing entangle- 
ment (^1 -> 0) and the loss of the visibility (i^ — > 0). 
Naturally, there is some optimally squeezed state, which 
give s the minimal value of A^9, as discussed in detail in 

We now calculate the sensitivity using the 

Gutzwiller approximation (l27l) for two wells, and check 
how it compares with the exact result (I37p . As shown in 
the Appendix [Cl in this case, both the one body Fisher 
information ([55]) and the coefficient C (IM)) can be eval- 
uated analytically and the outcome is 



m-N-A 6 



A'9 = 



1 



mN 



(38) 



where the tilde sign denotes the phase-squeezing param- 
eter and the visibility calculated with a product state. 
Although Equations (|57)) and (pS)) have identical form, 
we still have to verify whether the two- well entanglement 
of the fixed- TV state can be reproduced with a product 
state by increasing the on-site fluctuations. We model 
the coefficients of Eq. (P7|) with a Gaussian 



(39) 



where (n) = for a two-mode state. The on-site atom 
number variance is simply A^n = cr^, thus a coherent su- 
perfluid state can be modeled by setting the fluctuations 
at the Poissonian level, i.e. a = ^/N. To check if Eq. (|38|) 
will drop below the SNL when the fluctuations are super- 
Poissonian, we proceed as follows. For various values of 
a > a/TV we calculate the phase-squeezing parameter 
and the visibility j> and plot the sensitivity p8|) as a func- 
tion of To compare with the fixed- case, we find the 
ground state of the Hamiltonian (|26|) for various values 
of the ratio < 0. For every state found in this way, 
we calculate the phase-squeezing parameter and the 
visibility v and on the same plot draw the sensitivity ([37]) 
as a function of 5^. Figure [2] shows that not only Equa- 



tions (|55|) and ([57)) have identical form but also the quan- 
tum correlations giving sub shot-noise sensitivity can be 
perfectly mimicked using a product state (|27p with large 




FIG. 2. The sensitivity m - N ■ of the parameter estima- 
tion from the fit to the density of the interference pattern with 
M — 2 and A'^ = 200 calculated using the state with fixed- 
(open blue circles) and non-fixed number of atoms (solid black 
line) as a function of the phase-squeezing parameter. The sen- 
sitivities shown are normalized to the SNL. The dashed red 
line shows the ultimate precision of the phase estimation given 
by the QFI from Eq. ([15}. The horizontal black dotted line 
shows the shot noise limit. 



B. Fit sensitivity with M lattice sites 

For M > 2, when the integrals ([23]), ([33 and ([55)) 
cannot be evaluated analytically, we calculate the sensi- 
tivity ()23p numerically. We fix the average total number 
of atoms to be (N) = 5 x 10"* and for a given M we cal- 
culate the sensitivity ([^5)) for different values of the on- 
site atom-number fluctuations using the Gaussian model 

([55)) with (n) — In order to compare the results for 
different M, we normalize the sensitivity by the factor 
■|(Af^ — 1), that is, by the sensitivity achieved with the 
best possible measurement ([TTI at the shot-noise level. In 
Fig.[5]we plot m{N) x ^{M^-l)A'^e for M = 2, 4, 6, 8, 10 
and 12. The outcome is - with a high level of accuracy - 
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independent of M. We thus conclude that the sensitivity 
of the fit has a form similar to (|17p . i.e. 

^ (40) 



A^0 = 



1 



1 



mC{{N)) (M2- 1)^ 



where C,{{N)) is some function of the average total num- 
ber of atoms, which becomes larger for states with in- 
creased on-site atom number fluctuations, C > {N). 

m-N-A^e 




FIG. 3. (color online) The sensitivity m ■ N ■ A^6 of the 
parameter estimation from the fit to the density of the inter- 
ference pattern for M = 2, 4, 6, 8, 10 and 12 (different color 
solid lines) with a fixed average number of atoms equal to 
(iV) = 5 X 10*. The sensitivities shown are normalized to the 
SNL and multiplied by the factor i(M^ - 1). The dashed 
red line shows the ultimate precision of the phase estimation 
given by the QFI from Eq. (|15|) . The horizontal black dotted 
line shows the shot noise limit. 



vice for the measurement of spatially dependent forces. 
For a potential gx^ with j > 0, we showed that, opti- 
mizing over all possible output measurements, the cou- 
pling constant g can be measured with an uncertainty Ag 
scaling like , for a large number of lattice sites AI, 
and independent of the quantum state employed. We 
then showed that this scaling with the number of lat- 
tice sites can still be reached with a particular output 
measurement which is commonly employed with ultra- 
cold atoms: a least square fit to the interference pattern 
resulting from the overlap o the freely expanding atoms 
released from the confining potential. We also derived 
the improvements coming from entangling the quantum 
state. 

The device here studied is an example of a matter- wave 
multi-path interferometer. While the sensitivity bounds 
for a general two-path interferometer are already largely 
studied, the same analysis based on rigorous estimation 
theory results has not yet been performed for multi-path 
devices. Moreover, the potential of the implementation 
of these devices with atoms seems very promising, due to 
the growing control and tunability of the available trap- 
ping and probing tools, especially for ultracold atoms. 

Finally, we underline that the precision of the time-of- 
flight interferometer, where the phase is estimated from 
the fit of the averaged density to the measured interfer- 
ence pattern cannot be surpassed by adding a hopping 
term to the generating Hamiltonian The presence of 
the hopping - which induces the Bloch oscillations - is 
an extra complication, which changes the optimal states 
depending on the ratio between the tunneling energy to 
the phase related energy e, but does not introduce any 
further benefit. Some of these conclusions can be found 
in [s^l and will be discussed in detail elsewhere. 



In Fig. ini as the phase-squeezing parameter drops, the 
discrepancy between the curves calculated for different 
M grows. This is a spurious result due to the Gutz wilier 
approximation (l27t for the state \tp). When the aver- 
age number of atoms in the system is fixed and M is 
changed, the average on-site occupation varies. This is 
not a problem for as long as the width a of the Gaussian 
(p9)) is small as compared to (n). However, we increase a 
to mimic the growing entanglement in the system. Once 
a is comparable to (n), the change in the on-site popu- 
lation for different M becomes relevant. Therefore, by 
changing M we are actually changing the nature of the 
quantum state and, in particular, we cannot assume that 
the squeezing parameter remains fixed. One could avoid 
this problem by - rather than fixing (N) - fixing the on- 
site population (n) and increase the number of wells so 
that {N) =M{n). 

V. CONCLUSIONS 

We have analyzed the performance of ultracold bosons 
trapped inside an optical lattice as an interferometric de- 
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Appendix A: Calculation of the QFI 

In order to calculate the Fg, we need to derive an 
expression for the variance of the operator h. Using the 
definition ©, we get 

(Al) 
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We now make of the symmetry between the wells. The 
average value (n^) does not depend on index k and thus 
can be denoted as (n^) . The same argument concerns the 
inter- well correlation {fik-^fik^), which will be denoted as 
{fiifij), where i,j are indeces of any two sites. The above 
expression thus boils down to 

(h^) = {n') Y,e + (nSj) ^1^2 = (A2) 

k fci#fc2 

y (M + 1){2M + + {n^nj)^{M^ - l)(3Af + 2). 

Our next step is to show that {fiifij) can be written as 
a function of {v?) only. To this end, since the state is 
symmetric, we can take i = 1 and j — M and obtain 



UM-i) (A3) 



M 



Now, since it doesn't matter which two wells we correlate, 
in the sum instead of ni we can use an average over M—1 
wells as follows 



Using the result for the average value of h, we finally 
obtain 



Fq = 4A2/i = 1 (A/ + 1) Af 2 A^n, (A8) 
3 



where A^fi = {tt?) — (n)^ is the on-site variance. 



Appendix B: Calculation of the sensitivity 

In this Appendix we give details of derivation of co- 
efficient G from Eq. (p4)) which is used to evaluate the 
sensitivity ([^ using the product state (P71) . In order to 
calculate C, one first has to find the second order corre- 
lation function, 

G(2)(a;,y|6l) = {x\e)¥ {y\e)^ {y\e)i! {x\e)) . (Bl) 



When the field operator, which is introduced in Eq. 
is put into this definition, we obtain 



2 ni + . . . + riM-i 
Af- 1 



(ni+...+nM-i). G^^Hx,y\e) = \i' {■^) f |^(^) 



(B2) 



(A5) 

Finally, we exploit the fact that the total number of 
atoms is fixed and thus ni + . . . + um-i = N — um 
and obtain 

{niUM) ^^^l^Cfi , (A6) 



M ^ " Af - 1 

n 



M 



Af - 1 



(A7) 



M 

E 



i4){x)(ki-k4.) i<l>(y)(k2-k3) -A - 



k\ ,k2^k3^k4 — l 



«fci«fe2"fe3afc4>,(B3) 



where cjj^x) = 2^x + 9. The above sum cannot be per- 
formed in a straightforward manner, since special care 
has to be taken to distinguish cases when the indices 
ki . . . k4 in the average value (a|,^ a^^ a^g a^^ ) are equal or 
different from each other. A careful classification of all 
the possible terms gives 



X {M{a^a^aa) + {f{x + y) - M){a^a^f + {n)'^{f{x - y) + Af^) (B4a) 

+2 [fix) + f{y) - 2M] {a^a^a){a) + {h){a)^ [8M - 2M^ + (Af - 4) (/(x) + /(y)) - 2f{x - y) + 2g{x, -2/)](B4b) 
+ (a)2(a2) [4A^ - 2/(x) - 2/(y) - 2f{x + y) + 2g{x, y)] (B4c) 
+(a)4 [f{x)f{y) -{M- A){f{x) + /(y)) + f{x + y) + }{x - y) ~ 2g{x, y) - 2g{x, -y) + - QM] ] . (B4d) 



The function f{x) is defined as 

cos [M(j){x)\ 



cos [4>{x)] 
while the g{x,y) function is 



(B5) 



sin 


' Mip{x)' 
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sin 


\M<l>{y) 
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sin 


\ M{<j>(x)+^{y))] 
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sin 


\<t>{x)] 
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sin 


[0(^)1 
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sin 


' 4>{x)+<l,(y)' 
2 





(B6) 



r 



The coefficient C is defined as a double integral of 
the above second order correlation function with the 
one-body probabilities calculated at positions x and y, 
namely 

fdx fdyG^^\x,y\e)de\og[p{x\0)] delog[p{y\0)]. 



(B7) 

The one-body density is defined in Eq. ([30)) . We use 
the fact that the interference pattern consists of many 
fringes, so that the envelope varies slowly as com- 
pared to the oscillatory functions present in f, g and the 



10 



density p. Therefore, the integral (jB7l) can be rewritten 
as 

X de\og[p{m] de\og[p{(t>'\6)], (B8b) 

where we changed the variables x ^ (j) and y ^ (j)' . It 
can be demonstrated that 



27r 



a^iog = 0, 



a) 



where v is the fringe visibility calculated with a product 
state, 

.= l(i.)^l^^!M±i^^l(a)^ (C2) 

where we used the symmetry between the two wells, and 
the last step was possible because the coefficients C„ are 
taken to be real. We insert (jCip into ((33)) and obtain an 
analytical outcome 



Fi = 1 - v/r~ 



(C3) 



therefore all the terms from Equations (jB4p which do not 
depend on both (jj and 0' are zero. Using a property of 
the functions / and g 



d4> r d(j)' 



/(0 + 0'Hlog[p(0|0)] dg\0g[p{cp'\e)] 



' ^ r ^/(0-0'Hlog[p((/.|0)] de\og[p{^'\e)] 

and the same with g{(l), (j)') and g{(j)^ —4>') we finally get 
C = [{dY - {n? + 2(rj)(a)' - 2(a2)(a)2] x (BlOa) 



(BlOc) 



27r 27r 

+ [2(a2)(a)2_iv(a)2] x 

d(t> r d^'p'{^)p'{4>') 



27r 277 p(</.) 
which coincides with Eq. (p4|). 



(BlOd) 



Appendix C: Sensitivity for two wells 

First, we calculate the one-body Fisher information 
defined in Eq. ([33|) . For M = 2 the / function from 
Eq. ([HI]) reads 



/((/)) = 1 + t'COS0, 



(Cl) 



The integrals Ii from Eq. psp and /2 from Eq. ([36l) can 
be calculated in a similar way, giving 



h = h = (-2 + + 2vr^ 



(C4) 

i/^ 

As a result the coefficient C from Eq. ([M]) for M = 2 

reads 



C= ^(-2 + ip2^2Vr 



(C5) 



Finally, we note that for a symmetric state (Jy) — 0, 
thus the variance of the Jy operator calculated using a 
product state reads 



/~2\2 

- a 



(C6) 



which allows to rewrite Eq. (jC5l) in a following way 



C = — (-2 + z>2 + 2Vl-i>2) 



4 ^ -^y 



(C7) 



Combining ((C3)) and ((C7| as in Eq. 

1 



we obtain 



(C8) 



where is the phase-squeezing parameter calculated 
with a product state. 
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